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ABSTRACT 

We present results from the analysis of a mosaic of thirteen XMM-Newton pointings covering 
the Virgo Cluster from its center northwards out to a radius r ~1.2 Mpc (~4.5 degrees), 
reaching the virial radius and beyond. This is the first time that the properties of a modestly 
sized (M v j r ~ 1.4 x 10 14 M Q , IcT ~ 2.3 keV), dynamically young cluster have been studied 
out to the virial radius. The density profile of the cluster can be described by a surprisingly 
shallow power-law n e oc r~P with index = 1 .21 ±0. 12. In the radial range of 0.3r v ; r < r < r v ; r , 
the best fit temperature drops by roughly 60 per cent. Within a radius r < 450 kpc, the entropy 
profile has an approximate power-law form K oc r 11 , as expected for gravitationally collapsed 
gas in hydrostatic equilibrium. Beyond r ~ 450 kpc, however, the temperature and metallicity 
drop abruptly, and the entropy profile becomes flatter, staying consistently below the expected 
value by a factor of 2-2.5. The most likely explanation for the unusually shallow density 
profile and the flattening of entropy at large radius is clumping in the ICM. Our data provide 
direct observational evidence that the ICM is enriched by metals all the way to r2oo to at least 
Z = 0.1 Solar. 
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1 INTRODUCTION 

Sitting at the top of the mass spectrum of the virialized objects 
in the Universe, galaxy clusters serve as unique probes of cosmo- 
logical parameters (e. g. |Schuecker et al.|2003| |Allen et al.|2008| 
|Vikhlinin et al.|2 009; Mantz et al. 2010} and provide an approx- 
imately fair sample of the matter content of the Universe ( |White| 
|et al.|1993[|Carlberg et al.|1996> . 

Over the past 20 years, progress in understanding the astro- 
physics of galaxy clusters and their use as probes for cosmology 
was driven by X-ray observations. However, those observations 
only targeted the inner part of their volumes extending at most out 
to rsotf] where the X-ray emission is brightest; extrapolations were 
used to infer their properties further out. Until very recently, the 
thermodynamic properties of the cluster outskirts were not directly 
observable. 

Early work on the outskirts of clusters focused on determin- 
ing the surface brightness profiles, using ROSAT data in particu- 
lar (Vikhlinin et al. 1999; Neumann 20051. This provided density 
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estimates, but essentially no temperature information. Reliable in- 
formation regarding the temperature of the intra-cluster medium 
(ICM) at the virial radius (which approximately corresponds to 
r 2 oo) has become available only recently. Properties of the ICM 
in the outskirts of a handful of clusters, including PKS 0745-191 
( |George et a"ij|2009"l>, A bell 2204 ( |Reiprich et "aT1[2009), Abell 
1795 ( |Bautz et al.||2009| i, Abell 1413 flHoshino et al.||2010T >, and 
Abell 1689 (Kawahar ada et al.|20l"0l >, could be studied thanks to 
the low instrumental background of the Suzaku satellite. The most 
detailed radial profiles of thermodynamic quantities measured all 
the way out to ?'2oo were recently obtained for the Perseus Cluster 
(Simionesc u et al.|2011} . Among other results, these observations 
provide evidence that the ICM at large radii is highly clumped. 

However, all of these clusters studied with Suzaku at large 
radii to date are hot, massive, relatively relaxed systems. In order 
to have a complete picture of how galaxy clusters evolve within the 
context of hierarchical structure formation, we also need to under- 
stand less massive, less relaxed systems. Such systems are also crit- 
ical to understand for cosmological studies: a large fraction of the 
integrated Sunyaev-Zerdovich (SZ) power spectrum is expected to 
come from modestly sized clusters and groups, which are much 
more common than their more massive counterparts. 

The nearest, brightest, modestly sized, dynamically young 
cluster and therefore the ideal object to study for these purposes 
is the Virgo Cluster. The cluster is at a distance of 16 Mpc and has 
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an average temperature of ~2 keV. The bulk of the X-ray emission 
from this relatively cool system arises in the soft X-ray band, where 
the effective area of current X-ray mirrors is large and the signal-to- 
detector background contrast is optimal. Therefore, using a careful 
background treatment, we can study the thermodynamic properties 
and the chemical abundances of the low surface brightness regions 
of Virgo out to T2oo, even with the XMM-Newton satellite. 

|B6hringer et al.| (1994) used ROSAT Position Sensitive Pro- 
portional Counter (PSPC) observations from the ROSAT All-Sky 
Survey to study the structure of the Virgo Cluster. They discovered 
that the X-ray emission in general traces the galaxy distribution and 
confirmed 3 sub-clumps, centered on the galaxies M87, M49, and 
M86 respectively (reported earlier by Binggeli et al. 1987 1993| >. 
Shibat a~et al.| l |2001) derived an extensive temperature map using 
hardness ratio values from ASCA observations covering an area of 
19 deg 2 . For a mean temperature of kT = 2.3 keV, the scaling re- 
lations by |Arnaud et al.||2005"| predict M 20 o = 1.4 x 10 14 M Q and 
r 2 oo = 1-08 Mpc, which is a projected radius of ~ 3.9° in the sky. 

In this work, we study the structure of the Virgo Cluster us- 
ing a mosaic of XMM-Newton pointings, which cover it from the 
center northwards out to r 2 oo- The pointings are shown in Fig. [T] 
as circles over-plotted on the adaptively smoothed ROSAT PSPC 
image. Sect. [2] discusses the data reduction and extraction of the 
data products. In Sect. [3] we present projected surface brightness, 
temperature and metallicity profiles and determine the deprojected 
profiles of temperature and density. Based on these values we cal- 
culate the pressure and entropy profiles of the ICM. We also discuss 
in detail various possible sources of systematic errors and their im- 
pact on our results. In Sect. [4] we discuss the implications of our 
results and provide a summary and conclusions in Sect. [5] 

All the errors are given at lcr (68%) confidence level unless 
otherwise stated. We assumed the distance of Virgo as that of the 
central galaxy M87 (16.1 Mpc, |Tonry et al.|200l) . At this distance, 
1' corresponds to ~ 4.65 kpc. 



280 kpc 
60 arcmin 



r200 




Figure 1. Adaptively smoothed ROSAT PSPC image of the Virgo Cluster in 
the 0.5-2.0 keV band with the XMM-Newton pointing positions over-plotted 
as 12' radius circles. The color table corresponds to a logarithmic brightness 
scale. The background pointing lies 1.5° to the north of the northernmost 
pointing shown. M87 is the brightest cluster galaxy and sits in the centre of 
the Virgo Cluster, M86 is in the centre of an X-ray bright a sub-clump. The 
large elliptical galaxy NGC4477 and the spiral galaxy M88 are likely asso- 
ciated with sub-clusters and are marked with red crosses. The approximate 
T2oo of the Virgo Cluster is shown with a dashed circle. 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 Observations 

A total of fourteen XMM-Newton observations were used, the de- 
tails of which are listed in Table [T] The observations were carried 
out in June and July 2002, with the exception of the southernmost 
ones. Thirteen of the pointings are partially overlapping and cover 
a stripe in the sky with dimensions roughly 30' (East- West) by 4.5° 
(North-South). The combined image of these pointings using the 
MOS detectors is shown in Fig. [2] The northernmost observation 
located ~ 5.5° from the cluster center (not shown in Fig.[T]and|2| 
lies beyond the virial radius of Virgo, contains no significant ICM 
emission, and was used to determine the local X-ray background. 

2.2 Data analysis 

The data analysis was carried out using the XMM-Newton Science 
Analysis System (SAS) version 9.0.0. The filtering and extraction 
of the MOS data products and the subtraction of their instrumen- 
tal background was performed using the XMM-Newton Extended 
Source Analysis Software (XMH-ESAS) and methods described in 
Kuntz & Snowden| (2008}; |Snowden et"al~| (2008) ; |Snowden &| 
Kuntz| ^20 10) . The filtering and subtraction of the instrumental 
background for the pn data was based on that performed by |Wernerj 
|etal]((2008). 



In brief, to remove the soft proton flares, we used the 
mos-filter procedure in XMM-ESAS ( |Snowden & Kuntz|[2fjT0) . 
The MOS detector background was modeled and subtracted using 
the filter-wheel-closed (FWC) data together with data from the un- 
exposed corners of the CCDs. The total count rates and hardness ra- 
tios are determined in the unexposed regions of the CCDs for each 
source observation and a data base of all archived observations is 
searched for FWC data with similar count rates and hardness ratios. 
FWC data with the most similar count rates and hardness ratios are 
then scaled by the ratio of the source observation spectra from the 
CCD corners to the FWC spectra from the corners. These scaled 
FWC data are then used to produce the background files for the ap- 
propriate regions. The background modeling is done for each MOS 
chip individually. 

To remove the soft proton flares from the EPIC/pn detectors, 
we extracted light curves in the hard energy band, 10-14 keV, and 
created count rate histograms which have a roughly Gaussian peak 
at the nominal count rate. We excluded the time periods when the 
count rate exceeded the mean by more than 3cr. We repeated this 
procedure in the soft energy, 0.3-1 keV, band. The spectral proper- 
ties of the EPIC/pn detector background are relatively stable and 
the chip-to-chip variations are not significant. For EPIC/pn, we 
therefore used a stacked FWC data set which we again scaled by 
the ratio of the number of the counts in the unexposed corners of the 
source observation and the FWC data. We also tried to normalize 
by the ratio of the observed and FWC count rates in the 10-14 keV 
band, where practically all counts are instrumental and the number 
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Figure 2. Mosaic XMM-Newton image in the 0.6 - 3.0 keV energy range. The background pointing lies 1 .5° to the north of the northernmost pointing shown. 
Only the central region within the radius of 12' from the telescope's aim-point was used for each instrument. The image has been exposure-corrected, smoothed 
with a gaussian with the width of 3 pixels, and the point sources have been removed. The left panel shows the southern part of the mosaic and the right panel 
shows the northern part. The white circle marks the overlap in both panels. The central part of the cluster is clearly dominated by the emission of M87, the 
cold front discussed in Simionescu et al. (JOIOi is visible at ~ 20' to the north of the center. 



© 0000 RAS, MNRAS 000, 000-000 



4 O. Urban et al. 



Table 1. Summary of the observations. Columns give the observation identifier, starting date of the observation, XMM-Newton revolution of the observation, 
total and cleaned exposure times for MOS1, MOS2 and pn instruments respectively, and the coordinates of the pointings. The pointings in the table are ordered 
from the south to the north. 



obs ID 


date 
DD/MM/YY 




Exp. time flat, 
MOS1 MOS2 


pn 


Clean time [ks 
MOS 1 MOS 2 


pn 


da r°i 


DEC [°] 


0200920101 


1 1/01/05 


932 


95.36 


95.37 


_W 


74.87 


75.65 


- 


187.71 


12.36 


0106060101 


12/07/01 


291 


9.28 


9.28 


5 


8.75 


9.04 


5 


187.71 


12.73 


0106060201 


08/07/02 


472 


8.36 


8.37 


5 


4.65 


4.51 


1.7 


187.71 


13.06 


0106060301 


04/07/02 


470 


8.89 


8.89 


5.52 


6.87 


7.3 


5.1 


187.71 


13.39 


0106060401 


04/02/02 


470 


11.33 


11.35 


7.64 


10.54 


10.68 


6.77 


187.71 


13.73 


0106060501 


06/07/02 


471 


17.38 


17.42 


13.87 


14.3 


14.7 


11.58 


187.71 


14.06 


0106060601 


08/07/02 


472 


14.37 


14.37 


11 


9.38 


9.81 


8.52 


187.71 


14.39 


0106060701 


05/07/02 


471 


13.85 


13.87 


10.5 


4.89 


6.09 





187.71 


14.73 


0106061401 (6) 


13/06/02 


460 


8.23 


8.26 


4.88 


7.84 


7.92 


4.88 


187.71 


15.06 


0106060901 


10/06/02 


458 


15.03 


15.03 


11 


11.28 


11.38 


7.8 


187.71 


15.39 


0106061001 


06/06/02 


456 


12.39 


12.18 


8.89 


8.58 


8.45 


5 


187.71 


15.73 


0106061101 


09/06/02 


458 


15.93 


15.95 


11.98 


8.31 


8.43 


4.96 


187.71 


16.06 


0106061201 


09/06/02 


458 


17.57 


17.65 


13.52 


14.19 


14.36 


10.22 


187.71 


16.39 


0106061301 (c) 


10/06/02 


458 


15.99 


16.03 


12 


14.7 


14.39 


11.33 


187.71 


17.89 



(n) data from pn were not used for the central pointing 
replacement for flared observation 0106060801 
background pointing (see text) 



of source counts is negligible. Both methods of scaling produced 
consistent results. 

Table[T]shows the observation ID, date, revolution number, ex- 
posure times before and after cleaning for each instrument, and the 
coordinates of each pointing. For both EPIC/MOS and EPIC/pn, we 
check for residual soft proton contamination by comparing the 10- 
12 keV (10-14 keV for EPIC/pn) count rate per unit area within 
the field of view and in the unexposed corners. Soft protons are 
focused by the mirrors of the telescopes and they do not affect the 
unexposed areas of the detectors. Therefore, this comparison, in the 
hard energy band where the number of source counts is negligible, 
is a good way to detect a significant residual soft proton contam- 
ination (see Leccardi & Molendi 20081. Two of the observations 
(obs. ID 0106060701 and 0106061101) were affected strongly by 
flares and were not used in the spectral analysis due to residual soft- 
proton contamination. For observation ID 0106060201, the useful 
pn exposure time after cleaning was only 1.7 ks and the comparison 
of the exposed and unexposed areas of the detector in the hard en- 
ergy band indicates residual soft proton contamination. Therefore, 
it too has not been used in the further analysis. 

Images were extracted in the 0.6 - 3.0 keV energy range. We 
stacked the images for all three instruments and used the SAS task 
ewavelet to detect point sources down to 5cr significance. The 
stacked images were visually inspected to correct for false detec- 
tions, sources with larger dimensions, and point sources that were 
missed by ewavelet. For the central pointing, we used the list of 
point sources from an earlier analysis by Simio nescu et al.| ( |2010| >. 
For a z = 0.98 background cluster located about 1.3° north of M87 
( |Fassbender et al.|2010) , we excluded a region within its rsoo- A 
total of 580 point sources were excluded from the subsequent anal- 
ysis. 

To minimize systematic uncertainties associated with mea- 
surements at large off-axis angles (e.g. effective area calibration 
uncertainties), we have only used the data collected from within 
circular extraction regions of radius 12', centered on the aim points 
of the individual telescopes for each pointing. 



2.3 Spectral analysis 

We extracted spectra from all clean observations north of the central 
pointing. Previous results from these data for the cluster's central 
region are discussed by e .g.|Belsole et al. j2001}; Molendi ( 2002[>; 



et al.| ( |2006) ; |Simionescu et al.| ( |2007||2008l|2010> 



Bohringer et~al~| pOOT); |Matsushita et al.| j2002| |2003| >; |Werner| 



The spectra were extracted from a set of annular regions cen- 
tered on the active nucleus of M87 {a = 12h 30m 49.4s, S = 
12° 23' 28"). The widths of the annuli were varied to maintain 
a nearly constant number of instrumental background subtracted 
counts (10000) in each region. The SAS task rmfgen was used to 
create response matrices (RMF). Auxiliary response files (ARF) 
were created for each region using the task arfgen. 

We used XSPEC 12.6.8 jArnaud||1996| > to model the spec- 
tra. Fits were carried out in the 0.4 - 7.0 keV energy range us- 
ing the combined data from all three detectors. We excluded the 
1 .2- 1.8 keV interval due to the instrumental Al and Si lines. For the 
diffuse X-ray emission, we used a single-temperature apec model 
(Smith et al. 2001), which describes an optically thin plasma in 
collisional ionization equilibrium. The redshift was fixed to that of 
M87 (z = 4.36 x 10 -3 ). For each pointing, we fixed the Galac- 
tic absorption column density to the value determined by the Lei- 
den/Argentine/Bonn Survey (Kalberla et al. 2005 j). The metallici- 
ties are reported relative to the Solar abundances given by |Grevesse| 
& Sauval ( 1998). We used the extended C-statistic (which allows 
for background subtraction) in all the spectral modeling. 



2.4 X-ray foreground and background 

We extracted the cosmic X-ray foreground/background (CXFB) 
spectrum from the northernmost pointing using the combined data 
from all three detectors, after removing the point sources. We as- 
sumed a sum of three components in the CXFB model - an ab- 
sorbed power-law from the unresolved point sources {De Luca 
|& Mo lendi 2004}, an ~0.2 keV emission from the Galactic halo 
(Kuntz & Snowden 20001 and the Local Hot Bubble emission mod- 
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Table 2. Fore- and background model parameters. Fluxes are given in the 2- 
10 keV interval. The second column lists the photon index of the absorbed 
power-law component and the temperatures of the two foreground thermal 
components. GH stands for Galactic Halo, LHB for Local Hot Bubble. 



component 


r/Temperature [keV] 


flux 


-1 —2 j -2 

erg s cm deg 




power-law 

apecoH 

apec L HB 


1.44 ±0.03 
191 +0.005 

0.08 1 ' 


11.90+° * x 10-" 
(6.32 ± 0.32) x 10-' 6 

Mr*™ x u- 21 



fixed value 



eled by 0.08 keV thermal plasma (LHE0 Sidher et al.|l996||Kuntz 
|& Snowden"l|2000fr . During the CXFB model fitting, the normal- 
izations of the three components, the power-law index, and the 
temperature of the Galactic halo gas were left to vary indepen- 
dently. The metallicities of the foreground thermal plasmas were 
fixed to the Solar value. The absorption column density was fixed 
at Nh = 2.13 x 10 20 cnr 2 jKalberla et al.|2005| l. The resulting fore- 
ground/background parameters are summarized in Table|2] The ob- 
served flux of the power-law component is lower than the value 
reported by De Luca & Molendi (2004), which is to be expected 
given that the point sources detected and resolved in our analysis 
were excluded. 

All spectra used this same CXFB model, with the normaliza- 
tions scaled to account for the areas of the extraction regions. 

The north-polar spur region, with its elevated soft foreground 
emission, is relatively strong throughout the eastern parts of the 
Virgo Cluster. However, our pointings are in a direction that avoids 
this emission. 



3 RESULTS 

3.1 X-ray surface brightness 

In Fig. [3] we present the background-subtracted surface brightness 
profile of the Virgo Cluster obtained from the combined, exposure- 
corrected MOS images in the 0.6 - 3.0 keV band. It was extracted 
from 199 annular regions stretching from the center out to ~ 255' 
(projected distance ~ 1.2 Mpc). In the cluster center, out to ~ 8.5' 
= 40 kpc, the regions are 0.1' wide; at larger intermediate radii their 
width is ~1.35', and beyond 700 kpc the binsize increases to 2.7'. 
The gaps in the profile indicate missing data due to flared observa- 
tions. The background was assumed constant over the whole clus- 
ter, with a value equal to the average surface brightness of the outer 
background pointing. 

The profile shows a prominent discontinuity associated with a 
cold front at 90 kpc. This feature has previously been discussed by 
|Simionescu et al.H201(ty . Outside the cold front, the surface bright- 
ness distribution can be described by a power-law model S x oc r~" 
with index a = 1.34 ± 0.01. We verified that there is no signif- 
icant difference between the slopes of surface brightness profiles 
obtained for different low-energy bands. 

In the faint outskirts, beyond r ~ 700 kpc, the surface bright- 
ness starts to fluctuate significantly and the profile becomes more 
uneven. We note that the individual surface brightness bins in these 



2 The LHB emission is produced within the Galactic absorption column 
and therefore this component is unabsorbed. 
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Figure 3. Background- subtracted surface brightness profile in the 0.6- 
3.0 keV band. The gaps in the profile indicate missing data due to flared ob- 
servations. The virial radius, which we define as r2oo> lies at about 1 .08 Mpc. 
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Figure 4. Projected temperature and metallicity profiles. The virial radius, 
which we define as rjjoOi li es at about 1.08 Mpc. 



outer regions contain at least ~500 counts. The fluctuations are 
therefore formally highly statistically significant. 



3.2 Spectral results 

The results of our spectral fits are shown in Table [3] and in Fig. [4] 
Because we fitted the data using C-statistics, we do not have have 
a goodness of fit indicator similar to the reduced x 1 - In order to 
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Figure 5. Examples of the analyzed spectra from ~0.5r2oo and from ~r2oo- Data from all three EPIC instruments - MOS 1 (black), MOS 2 (red) and pn (green) 
- were fitted simultaneously. Left: Data with over-plotted CXFB model. Bottom part of each plot shows the ratio of the total to CXFB flux, thus indicating the 
amount of cluster emission. Right: Spectra with over-plotted best fit model, which includes the cluster emission. 



confirm that the fits are acceptable, we re-binned the data to a min- 
imum of 50 counts per bin and evaluated our best fit model using 
X 2 statistics. All the reduced^ 2 values for our best fit models are in 
the range of 0.86-1.38. These values show that, given the present 
systematics, our fits are acceptable. The goodness of fit estimated 
from the;^ 2 statistics does not show a systematic trend with radius: 
the fits at large radii are, on average, as good as the fits at small 
radii. 

Both the temperature and metallicity profiles show discon- 
tinuities at the cold front (r ~ 90 kpc), previously studied by 
|Simion escu et al. (2010) and modeled by Roediger e t al.| < [20To) >. 
This feature can be explained by the sloshing of the dense gas as- 
sociated with the cooling core in the gravitational potential of the 
cluster. 

Beyond a radius r ~ 450 kpc, the temperature and metal- 
licity both drop abruptly. At the virial radius, the measured tem- 
perature reaches kT = 1.08^JJj keV and the metallicity is Z = 



0.11 ± 0.02 Solar. This is the first time that the metallicity at r2oo 
has been measured with high (>5.5cr) statistical significance. This 
is possible in Virgo due to its proximity and the low ICM tempera- 
ture, for which the Fe-L line emission is strong. At ~1 keV, where 
the Fe-L emission peaks, the effective area of the XMM-Newton 
mirrors is large and the signal-to-detector background contrast is 
optimal. 

The strength of the Fe-L signal can be seen in Fig. [5] where we 
show the spectra from the radius 522-580 kpc or r ~ 0.5 r 2 oo and 
for our outermost source extraction region at 1029-1198 kpc, or 
r ~ r 20o- The solid lines on the left show the background models for 
the given spectra; the bottom panels show the ratio of the measured 
X-ray flux to the CXFB flux. We can see clear excess signal above 
the background, representing the cluster emission. A large part of 
the detected signal from near the cluster's virial radius arises in the 
0.8-1.2 keV range, and is associated with Fe-L line emission. The 
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Table 3. Projected temperatures, metallicities, and normalizations obtained 
from spectral fitting. The first and the second columns give the inner and 



r (Mpc) 

Figure 6. The deprojected electron density («_), entropy (K), and pressure 
(P) profiles. For the density profile, we overlay the best fit power-law model, 
« c oc r~' 2 . On the entropy profile we over-plot K oc r 1 1 , which is expected 
for gravitationally collapsed gas in hydrostatic equilibrium. On the pressure 
profile, we over-plot the extrapolated average profile determined for a sam- 
ple of clusters within ~0.5r2oo by |Arnaud et al.| (2010i. The virial radius, 
which we define as r2oo> lies at about 1.08 Mpc. 



right panels show the same data with the over-plotted best fit model 
including the cluster emission. 

Deprojected profiles of the temperature and density were de- 
rived from the combined set of spectra using the pro jet model in 
XSPEC. In order to ensure the stability of the fit, we linked some 
of the parameters in neighboring annular regions. The deprojected 
density profile is shown in the top panel of Fig. [6] The measured 
electron density approximately follows a power-law model n c oc r~& 
with index p= 1.21 ±0.12. 

Based on the deprojected values of the temperature and 

density, we also calculated the entropy = and pressure 

(p = n e kT) profiles. These are also shown in Fig. [5] On the entropy 
profile we over-plot K oc r lA , which is expected for gravitationally 
collapsed gas in hydrostatic equilibrium (Tozzi & Norman 2001 
|Voit|2 005 ). Out to r ~ 450 kpc, the entropy profile follows the ex- 
pected power-law shape, but further out it becomes flatter, staying 
consistently below the expected value by a factor of 2-2.5. On the 
pressure profile, we over-plot the average profile for a sample of 
clusters observed within ~0.5r 2 oo ^Arnaud et al.| 2010l extrapolated 
to larger radii. Our measured profile shows clear departures from 
this average shape: a deficit is observed at ~200 kpc, and enhance- 
ments are seen at ~400 kpc and ~600 kpc. Because the Virgo Clus- 
ter is relatively unrelaxed, the results of the deprojection analysis, 
which assumes spherical symmetry, should be treated with caution. 
The systematic uncertainties due to the departures from spherical 
symmetry are significantly larger than the measurement errors. 

The last data point in the deprojected electron density profile is 
most likely overestimated due to unaccounted emission from out- 
side of this radius, biasing the pressure in the last radial bin high 
and the entropy low. A small ringing effect from this artifact will 
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also be projected onto the next bin inwards. Such artifacts do not 
affect our main results on the density and entropy profiles. The ob- 
served pressure enhancement at 400 kpc coincides with a region of 
increased X-ray surface brightness, a possible sub-group surround- 
ing the massive elliptical galaxy NGC4477, lying just outside of our 
XMM-Newton extraction region (indicated in Fig. [TJ. The pressure 
increase at 600 kpc coincides with an X-ray bright region surround- 
ing the massive spiral galaxy M88 (see Fig. [TJ. Projection effects 
associated with the gas of these sub-groups could be responsible 
for the observed features in the pressure profile. 

The flattening of the entropy profile is robustly measured and 
particularly interesting. The observed entropy beyond the radius of 
~450 kpc is consistently a factor of 2-2.5 lower than the expected 
value, implying that the observed density is higher and/or the ob- 
served temperature is lower than expected. 



3.3 Systematic errors 

In order to test the robustness of our results, we fitted the data ob- 
tained by the three XMM-Newton EPIC detectors separately. The 
temperature results are shown in the left panel of Fig. [7] While the 
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Figure 7. Left panel: Temperature profiles for the separate EPIC instruments - MOS 1 (black), MOS 2 (red) and pn (blue). EPIC/pn was not used in the 
analysis at the radius at which the pn value is shown with dashed lines. Right panel: Ratio of the flux of the ICM versus the flux of the CXFB for the annular 
regions. Fluxes were determined in the 0.4 - 3 keV range. 
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Figure 8. Projected temperatures (upper panel) and normalizations (lower 
panel) with over-plotted profiles obtained by varying CXB parameters (see 
the text for details). We show the bracketing values due to the uncertainties 
on the normalization of the power-law background component (norm pow ), 
normalization of the Galactic halo emission (normoHX normalization of the 
LHB emission component (nonriLHB): and the temperature of the Galactic 
halo emission component (temp GH ). 



results obtained by the MOS 1 and MOS2 detectors agree remark- 
ably well at almost all radii, the temperature values obtained by 
EPIC/pn are systematically higher. The largest inconsistency be- 



Figure 9. Projected temperatures (upper panel) and normalizations (lower 
panel) with over-plotted profiles obtained by varying the instrumental back- 
ground (see the text for details). 



tween the MOS and pn detectors is in the one but last annulus 
(shown with dashed lines in Fig.|7J. We can find no reason for this 
inconsistency - it appears not to be due to residual soft proton con- 
tamination or technical issues. On the grounds of their gross incon- 
sistency with the two independent MOS detectors (which are fully 
consistent with each other), at this radius we have excluded the 
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EPIC/pn data from the analysis. With the EPIC/pn data included, 
the best fit temperature in this region is 1.46 ± 0.17 keV. 

However, we emphasize that the results are in a relatively good 
agreement at most radii and the general trend with the temperature 
decreasing at large radii is robust. 

3.3.1 Cosmic X-ray fore- and background 

The right panel of Fig. [7] shows the radial dependence of the ratio 
of the ICM/CXFB fluxes in the 0.4 - 3.0 keV range. The ratio drops 
below unity beyond ~ 800 kpc. However, as we show in the left- 
hand side of Fig. [5] the signal at ~1 keV is relatively strong even in 
the outermost annulus. 

We determine the variation of the soft foreground components 
within the area covered by our pointings using the images obtained 
in the ROSAT All Sky Survey (RASS) in the soft 0.1-0.3 keV 
(R1+R2) band. Most of the LHB and GH foreground is emitted 
in this energy band. Studying the point-to-point variation in these 
ROSAT all-sky maps we determined a conservative upper limit of 
30 per cent for the variation of the LHB and Galactic halo fluxes 
along our investigated region. The Galactic halo temperature was 
found using the ROSAT hardness ratio maps to vary by up to 20 
per cent of its value. Due to the exclusion of point sources down 
to a faint level, the root mean square variation of the normaliza- 
tion of the power-law component is expected to be smaller than the 
statistical uncertainty with which its value has been determined in 
the background pointing. Thus we varied this parameter by its 3cr 
statistical error. 

These uncertainties were used to determine the bracketing val- 
ues for the temperature and normalization profiles shown in Fig. [8] 
The background parameters were up- and downscaled one at a time, 
leaving the others with their original values, and the model was re- 
fitted. Although at radii larger than ~500 kpc the systematic errors 
dominate over the statistical uncertainties, the values remain well 
determined and the trends in the profiles are robust. 

3.3.2 Instrumental background 

The uncertainty in the normalization of the instrumental back- 
ground is 3-5 per cent. These uncertainties arise mainly from the 
Poisson error associated with the scaling factor which is used to 
account for the variation of the instrumental background level with 
time. The scaling factors of all three detectors were simultaneously 
raised or lowered by lcr to achieve a conservative estimate. The 
bracketing values for the temperature and normalization profiles 
obtained in this fashion are shown in Fig. [9] Again, the trends sug- 
gest that the results are robust. 



4 DISCUSSION 

We have presented the first well-constrained measurements of the 
ICM properties for a modestly sized, dynamically young galaxy 
cluster out to large radii. Here we compare our findings with the 
results obtained for hotter, more massive clusters, and with pre- 
dictions based on numerical simulations. Finally, we discuss the 
implications. 



4.1 The unrelaxed outskirts of a forming cluster 

The surface brightness profile along the northern direction of the 
Virgo Cluster is unusually shallow. For 100 kpc < r < 1200 kpc it 



follows a power-law shape Sx K r~ a with index a = 1.34 ± 0.01, 
which is significantly smaller than a ~ 3 for Abell 1795 (Bautz| 
|et al.|20 09) or a = 3.6 for a sample of massive luminous clusters 
obser ved with Chandra between 0.4—0.7 r2oo, (Ettori & Balestra 



2009 1, Simulations of low mass clusters (M vir < 10 is M G ) by 



Ron- 



carel i et al. (2006) predict the ratio between the surface bright- 
nesses at j-200 and 0.3r2oo to be ~ 8 x 10~ 3 , which is more than 
an order of magnitude smaller than the value of 0.16 measured in 
Virgo. The emission in the outskirts of this young forming cluster 
is therefore considerably brighter than expected. 

The derived density profile follows a power-law shape n e oc 
r~P with index = 1.21 ± 0.12. For comparison, the density profile 
of the Perseus cluster has /? = 1.68 ± 0.04 (Simionescu et al.|201 1| > 
and that of Abell 1795 is even steeper (/3 = 2.27±0.07). Simulations 
by |Roncarelli et a l. (2006) predict that the gas density in the cluster 
outskirts (out to ~ 1.2 r vir ) can be described by a power-law model 
with an index [} = 2.4 - 2.5. 

We measure a ~60 per cent drop in temperature between 
0.3 - lr 2 oo, which is consistent with previous observations of some 
more massive relaxed clusters with Suzaku (Reiprich et al. 2009; 



|Bautz et al.||2009"l |Hoshino et aTJ|20l0l |Simionescu et al.||2011| >. 
In PKS 0745-191 the temperature drops by a factor of 4 ( |George| 
|et al.|2 009). Simulations (Ronca relli et al.|2 006i predict a drop of 
~40 per cent, which is also consistent with the analytic result by 
IQstriker et aL1 < [2005> . 

As we show in the middle panel of Fig. [6] in Sect. |3.2| out- 
side r ~ 450 kpc the measured entropy is consistently lower, by 
a factor of 2-2.5, than values expected from hydrodynamic simu- 
lations of gravitational ly collapsed gas in hydrosta tic equilibrium, 
which predict K oc r 1 1 (Tozzi & Norman|200i]|Voit|2005). Flatten- 



ing entropy profiles have also been seen in PKS 0745-191 i George 



ct al 



2009 1 and especially in the Perseus Cluster ( Simionescu et al 



201 1 1; however, in these hotter systems the entropy profile starts to 
deviate from the expected shape at larger scaled radii. 

|Simionescu et al.| ( |201 1| > show that the break in the measured 
entropy profile, and a simultaneous increase in the apparent gas 
mass fraction above the cosmic mean in the outskirts of the Perseus 
Cluster, are most likely caused by clumping in the ICM. The quan- 
tity that we infer from the observed deprojected emission measure 
is the average of the square of the electron density (« 2 ), not the av- 
erage of the electron density squared (n) 2 . Therefore, if the gas is 
clumpy, its non-uniform density will lead to an overestimate of the 
average electron density calculated from the emission measure. If 
the dense clumps are in pressure equilibrium with the surrounding 
gas, then they will be cooler than the ambient medium. However, 
it is more likely that they will be falling into the cluster, moving 
through the ICM and therefore partly confined by ram-pressure. 
Depending on the ratio of ram-pressure to thermal pressure support 
of the clumps, the average temperature in the multiphase regions 
will be underestimated because denser and cooler blobs have larger 
volume emissivity. Therefore, increasing clumping as a function of 
radius will result in a flattening of the observed surface brightness 
and density profiles and a steepening of the observed temperature 
profile. The over-estimated density and underestimated tempera- 
ture at large radii will then combine into a flattening of the entropy 
profile, as observed. 

In order to test further whether clumping can be responsible 
for the observed features in our profiles, we simulated an XMM- 
Newton spectrum of a two phase ICM, where both phases have a 
metallicity of Z = 0.3 Solar. A hotter 2 keV component is used to 
represent the ambient ICM and a 1 keV phase represents the clumps 
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of cooler plasma. The emission measures were adjusted so that, as- 
suming thermal pressure equilibrium, the cooler clumps will have a 
volume filling fraction of / = 0.2 (additional ram-pressure support 
will make the volume filling fraction of the clumps smaller). The 
simulated spectrum contains all the observed background compo- 
nents in the Virgo Cluster and has a similar statistical quality to our 
observed spectra from large radii. When fitting this spectrum with 
a single temperature thermal model we obtain a good fit, with a 
temperature of kT = 1.24 keV and a metallicity of Z = 0.16 So- 
lar, underestimating the true metallicity by about a factor of 2. We 
also overestimate the density of the ICM by a factor of 1.6. The 
entropy inferred from this fit is about 2.2 times lower than the true 
entropy of the ICM, indicating that such multi-phase medium could 
indeed be responsible for the flattening entropy profile. The pres- 
ence of cooler, denser blobs also explains the abrupt decrease of 
the temperature and observed metallicity beyond 450 kpc, where 
the entropy starts to deviate from the expected power-law shape. 

Clumping therefore appears to provide a natural explanation 
for the properties of the observed density, temperature, entropy, 
pressure, and metallicity profiles in the Virgo Cluster. We note 
that the real ICM is, however, probably multi-phase with different 
clumps having slightly different temperatures and the clumps being 
partly ram-pressure confined. 

While X-ray observations provide the measurement of (« 2 ), 
SZ observations depend linearly on n. Therefore, in regions with 
clumped ICM, we expect to see a discrepancy between the X-ray 
and the SZ signals. If there is no bias in measuring the temperature 
(e. g. the clumps are completely ram-pressure supported and have 
the same temperature as the ambient plasma) then this difference 
will provide an independent measurement of the clumping factor. 
Temperature biases arising from the multi-phase structure of the 
ICM may complicate the interpretation. 

4.2 Metals at large radii 

We present a metallicity profile spanning from the cool core region 
all the way out to r 2 oo- Our metallicity measurement is dominated 
by the line emission of the Fe-L complex and therefore primar- 
ily represents the Fe abundance. Our data provide direct observa- 
tional evidence that the ICM is enriched by metals all the way to 
rooo- In the range of -0.5-1 r yK , the observed metallicity profile is 
consistent with being flat at a value of about Z = 0.1 Solar. Pre- 
viously, the ICM metallicity at large radii has been obtained by 
|Fujita et al.H2008) from a large region spanning 0.5-1 /"200 in the 
compressed area between the clusters Abell 399/401 at the onset of 
a merger. They found that the ICM at large radii has been enriched 
to Z ~ 0.3 Solar (with respect to the Solar abundanc es of | Grevesse| 
|& Sauval|1998} . More recently, |Simionescu et al.| ( |2011) showed 
that the outskirts of the Perseus Cluster have also been enriched to 
a similar metallicity of Z ~ 0.3 Solar. 

The gas within our analyzed region is, however, most probably 
multi-phase, consisting of multiple temperature components. It has 
been shown that the metallicity is often sensitive to the modeling 
of the underlying temperature structure and, if multi-temperature 
plasma is modeled with a single-temperature spectral model, the 
metallicity of gas can be significantly underestimated (e.g. |Buote| 
|2000| see also our simulation in Sect. |4.1[ >. This is especially true 
around the temperature of 1 keV, where the Fe-L complex is very 
sensitive to the underlying temperature structure. Our measured 
metallicity could thus be an underestimate of the true metallicity 
and should be considered a lower limit. 

Nevertheless, using the best fit metallicity values and assum- 



ing that both the gas density and the metal distribution are homo- 
geneous, we calculate the cumulative Fe mass of the Virgo clus- 
ter. The total Fe mass obtained this way is ~ 4 x 10 9 M a within 
0.1 - 1^200, approximately half of which resides outside of 0.5i?200- 
Because the metallicity is likely biased low and neither the gas den- 
sity nor the distribution of metals are likely to be homogeneous, this 
metal mass is almost certainly an underestimate. We conclude that 
the total mass of Fe beyond 0.5r200 is at least 2 x 10 5 M Q . 

The two dominant mechanisms by which metals pollute the 
ICM are galactic winds and ram-pressure stripping. Galactic winds, 
driven by the thermal energy of a large number of supernova ex- 
plosions (e.g. Heckman 2003), are expected to dominate the en- 
richment of the intergalactic medium at redshifts of z = 2-3, when 
the bulk of the star formation happened. Part of this pre-enriched 
inter-galactic medium subsequently fell into the gravitational po- 
tential of clusters where it got shock heated and today constitutes 
the ICM. 



|Fujita et al. (2008) interpret the high metallicity in the out- 
skirts of Abell 399/401 as evidence for early enrichment by galac- 
tic super-winds. They argue that in order to exceed the pressure 
of ~ 2 x 10~ 12 Pa necessary to start stripping the galaxies |Fu-| 
[jita & Nagashima 1999} their velocities would have to exceed 
~ 2000 km s -1 . Following this argument, the relative velocities be- 
tween the galaxies of the Virgo Cluster and the ICM would have 
to exceed the unrealistically high value of ~ 3000 km s~'. How- 
ever, Hi imaging surveys have revealed several galaxies falling 
into the Virgo Cluster that display long H 1 tails at projected radii 
of 0.6-1.2 Mpc from the center of the cluster (e.g. |Kenney et al.| 
|2004l|Chung et al.|2007[|Kantharia et al.|2008| >. There are further- 
more several examples which suggest that stripping also operates 
in lower density, lower velocity dispersion environments of poor 
clusters and groups of galaxies (Kant haria et al.|2005||Levy et al.| 
|2007| >. Recent simulations (e.g. |Roediger & Bruggen|2007||2008| > 
show that galaxies can get stripped already at large radii with low 
ICM densities due to continuous or turbulent/viscous stripping (see 
|Nulsen|[l982) . In principle, all the metals outside 0.5^200 could 
be supplied by the complete stripping of ~150 galaxies, assum- 
ing 10 10 M Q of stripped gas, enriched to the Solar metallicity, per 
average galaxy. Stripping would, however, produce a radially de- 
creasing metallicity profile at the large radii, while the observed 
profiles in both Virgo and Perseus are consistent with being flat. 



The outskirts of Virgo have three times lower apparent metal- 
licity than the outskirts of the Perseus Cluster and the compressed 
region between the clusters Abell 399/401. A similar ultra-low 
metallicity of Z = 0.15 Solar has only been measured in the group 
of galaxies NGC 5044 at radii 0.2-0.4r vir JBuote & Lewis|2004| >. 
Given that the stellar over ICM baryon fraction of clusters is de- 
creasing with the increasing cluster mass (e.g. Gonzalez et al.|2007| > 
explaining the higher metal content of the more massive clusters 
would be challenging. However, this difference may be primarily 
due to the multiphase gas biasing the Fe abundance low in cooler 
systems more strongly than in hotter systems. While the metallicity 
of Virgo at large radii has been determined based on the Fe-L lines, 
the metallicities of hotter clusters, kT > 3.5 keV, are determined 
based on their Fe-K lines. The metallicity of the hotter ICM could 
be biased both toward higher and lower values, depending on the 
temperature structure. The bias is however much smaller than in the 
cooler systems, only of the order of up to ~30 per cent (e.g. |Ras!a| 
|et al.|2008|[Simionescu et al.|2009||Gastaldello et al.|2010| >. 
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5 CONCLUSIONS 

We present the results of the analysis of a mosaic of XMM-Newton 
observations of the modestly sized, dynamically young Virgo Clus- 
ter from the center to the north out to ~ 1.2 Mpc. We detect X-ray 
emission all the way out to this radius, which allows us to study 
the ICM properties in unprecedented detail. The surface brightness 
and density profiles are significantly shallower than predicted by 
simulations and measured in other more massive systems. The tem- 
perature is measured to drop by ~60 per cent between 0.3 - 1 r vil . 
Beyond the radius of ~450 kpc the temperature and metallicity drop 
abruptly and the entropy profile deviates from the power-law shape 
K oc r [A expected for gravitationally collapsed gas in hydrostatic 
equilibrium l |Tozzi & Norman|200l[|Voit|2005^ At these radii, the 
entropy is consistently lower than expected by a factor of 2-2.5. 

The most likely explanation for the unusually shallow den- 
sity profile, the flattening of the entropy profile, and sharp drops in 
temperature and metallicity at r > 450 kpc is the onset of signif- 
icant gas clumping, which then increases as a function of radius. 
The clumping in the Virgo Cluster becomes significant at smaller 
scaled radii than has been detected in the more massive Perseus 
Cluster (Simionescu et al.|20lT) . 

Our data provide direct observational evidence that the ICM is 
significantly enriched with metals all the way to r 2 oo- The measured 
metallicity profile flattens out in the cluster outskirts and reaches 
Z = 0.1 1 ± 0.02 Solar at r^oo- The metallicity at large radii is, how- 
ever, likely underestimated significantly because of the clumping 
and therefore this measured value should be considered a lower 
limit. 

The Virgo Cluster, being the closest galaxy cluster to us, pro- 
vides a unique opportunity to study large-scale structure formation 
in detail as it happens. The similarities and differences between 
the results for the Virgo Cluster presented here and for the Perseus 
Cluster discussed by Simionescu et al. (201 1 1 point to a wealth of 
interesting physics at large radii of clusters that is just beginning 
to be probed. Understanding this physics presents an exciting chal- 
lenge for numerical simulations. 
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